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ABSTRACT 

A  class  of  approximations  {Ss,m}  ^  periodic  function  /  which  uses  the  ideas  of  Fade,  or 
rational  function,  approximations  based  on  the  Fourier  series  representation  of  /,  rather  than 
on  the  Taylor  series  representation  of  /,  is  introduced  and  studied.  Each  approximation  Sn,m 
is  the  quotient  of  a  trigonometric  polynomial  of  degree  N  and  a  trigonometric  polynomial 
of  degree  M.  The  coefficients  in  these  polynomials  are  determined  by  requiring  that  an 
appropriate  number  of  the  Fourier  coefficients  of  Snm  agree  with  those  of  /.  Explicit 
expressions  are  derived  for  these  coefficients  in  terms  of  the  Fourier  coefficients  of  /.  It  is 
proven  that  these  “Fourier- Fade”  approximations  converge  point-wise  to  +  f{x~))/2 

more  rapidly  (in  some  cases  by  a  factor  of  1/Ar*^)  than  the  Fourier  series  partial  sums 
on  which  they  are  based.  The  approximations  are  illustrated  by  several  examples  and  an 
application  to  the  solution  of  an  initial,  boundary  value  problem  for  the  simple  heat  equation 
is  presented. 
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tract  No.  NASl-19480  while  the  author  was  in  residence  at  the  Institute  for  Computer  Applications  in 
Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681. 
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1 .  Introduction.  Fourier  series  are  used  widely  in  many  branches  of  applied  math¬ 
ematics.  For  example,  they  are  often  used  together  with  separation  of  variables  to  con¬ 
struct  analytical  solutions  to  boundary  value  problems  for  differential  equations  and 
with  a  variety  of  spectral  methods  to  find  approximate  solutions  to  these  problems  nu¬ 
merically.  For  practical  purposes,  approximate  solutions  to  these  problems  are  often 
obtained  using  only  a  finite  number  of  the  terms  in  a  Fourier  series.  This  trunca¬ 
tion  procedure  often  leads  to  nonuniformly  valid  approximations.  For  example,  when 
the  function  being  approximated  has  a  point  of  discontinuity,  the  Gibbs  phenomena  is 
present.  The  “oscillations”  caused  by  this  phenomena  typically  propagate  into  regions 
away  from  the  singularity,  and,  hence,  degrade  the  quality  of  the  partial  sum  approx¬ 
imation  in  these  regions.  Even  if  the  function  being  modeled  is  analytic,  but  has  a 
region  of  large  slop>e  or  curvature,  there  can  be  significant  oscillations  in  the  partial 
sums  outside  these  regions,  which  again  degrade  the  quality  of  the  approximation. 

Recently,  Gottlieb  and  Shu  [6]  and  Gottlieb,  et.al.  [7]  have  proposed  a  way  of 
overcoming  the  Gibbs  phenomena.  Their  technique  involves  the  construction  of  a  new 
series  using  the  Gegenbauer  polynomials  C^{x).  For  a  function  /  that  is  analytic  on 
the  interval  [—1, 1],  but  is  not  periodic,  they  prove  that  their  technique  leads  to  a  series 
which  converges  exponentially  to  /  in  the  maximum  norm.  To  do  this,  they  require 
that  the  parameter  A,  which  appears  in  the  weight  factor  (1  —  grows  with 

the  number  of  Fourier  modes  considered.  As  we  shall  demonstrate  below,  the  family 
of  approximations  we  shall  introduce  can  be  defined  explicitly  in  terms  of  the  known 
Fourier  coefficients.  This  will  prove  to  be  particularly  useful  for  certain  applications  (as 
we  shall  demonstrate)  when  the  Fourier  coefficients  are  themselves  functions  of  one  or 
more  other  variables.  Although  the  approximations  we  shall  define  do  not  “eliminate” 
the  Gibbs  phenomena,  they  do  mitigate  its  effect,  as  we  shall  show.  This  is  especially 
true  outside  a  “small”  neighborhood  of  a  point  of  discontinuity  of  /,  where,  for  practical 
purposes,  the  “unwanted”  oscillations  can  essentially  be  eliminated. 

To  fix  notation,  we  let  /  (x)  be  a  piece-wise  smooth,  27r-periodic  function.  Then 
we  can  Msociate  with  /(x)  its  Fourier  series  5(x)  defined  by 

OO 

(1)  5(x)  =  ao/2  -f  cos(nx)  -t-  K  sin(nx)), 

n=l 

(2)  an  =  {l/Tr)J  /(x)  cos(nx)dx,  b„  =  {l/ir)J  /(x)  sin(nx)dx, 

n  =  0, 1, 2, ... .  (For  some  of  our  formulas  below,  it  is  convenient  to  regard  the  coefficients 
{a„,  bn}  as  being  defined  by  (2)  for  negative  as  well  as  positive  values  of  n.  Thus,  a_„  = 
On  and  6_„  =  — 6„  for  all  integers  n.)  It  is  well  known  (see  [2],  for  example)  that  5(x) 
converges  to  /(x)  at  each  point  x  where  /  is  continuous  and  to  (/(x‘'')-t-/(x“  ))/2  at  each 
point  X  where  /  is  not  continuous.  For  practical  purposes,  /(x)  is  often  approximated 
by  its  Fourier  partial  sum  Snjo{x)  defined  by 

N 

(3)  SNfi{x)  =  ao/2  -f  53(a„  cos(nx)  -|-  bn  sin(nx)). 

n=l 
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We  now  define  a  class  of  approximations  to  /  which  uses  the  ideas  of  P2wile  approx- 
imants  (see,  e.g.,  [1]),  except  that  the  approximations  are  based  on  the  Fourier  series 
representation  of  /,  rather  than  on  the  Taylor  series  representation  of  /.  Thus,  for  any 
two  non-negative  integers  N  and  Af,  we  define  a  family  of  “Fourier- Fade”  approxima¬ 
tions  5iv,A/(ar)  by 

fA\  q  +  E!)Li(AnCos(nx)  -|-  Bn  sin(ng)) 

1  +  Em=i  (C'm  cos(ma;)  -t-  £>„  8in(mx))  ’ 

The  2Af  -I-  2N  -f  1  constants  { A„,  B„,  Cm,  Dm)  which  appear  in  the  definition  of  Sjv,m 
are  determined  by  the  condition  that  2Af  ■+■  2N  -t- 1  of  the  Fourier  coefficients  of  Sf/,M 
agree  with  those  of  /,  i.e. 


(5) 

(iM 

j  5A/,Af(a:)cos(jx)di  =  aj,  0<j  <J 

(6) 

(I/tt) 

w 

f  sin(A:x)dx  =  bk,  1  <  k  <  K 

where  J  +  K  = 

2Af  -f-  2N. 

The  idea  of  constructing  Fade  approximations  based  on  series  representations  of 
functions  other  than  the  classical  power  series  representation  has  been  suggested  and 
studied  by  several  other  investigators.  For  example,  Maehly  [9]  has  suggested  an  ap¬ 
proach  to  determine  the  coefficients  in  rational  approximations  based  on  Chebyshev 
series,  and  this  approach  has  since  been  described  in  numerical  terms,  with  examples, 
by  Ralston[10]  and  by  Fike[5].  An  account  of  the  basic  theory  of  such  approximations 
has  been  given  by  Cheney  [3],  who  considered  the  more  general  case  of  expansions  in 
terms  of  a  basis  whose  elements  satisfy  relations  of  the  form  ipupj  = 

Clenshaw  and  Lord  [4]  have  reviewed  rational  approximations  based  on  Chebyshev  series 
and  present  a  well-conditioned  method  for  constructing  the  elements  of  a  “Chebyshev- 
Fade”  table.  The  general  form  of  a  Fourier- Fade  representation,  such  as  equation  (4), 
has  also  been  suggested  by  Cheney.  He  has  proven  the  existence  of  a  rational  trigono¬ 
metric  function  which  “best”  approximates  a  continuous  function  /,  but  does  not  discuss 
any  detailed  algorithms  for  the  computation  of  the  coefficients,  rates  of  convergence  of 
the  approximations,  etc. 

In  the  sections  2-4  below  we  consider  the  special  case  of  equation  (4)  when  M  —  1. 
In  particular,  in  sections  2-3  we  develop  explicit  formulas  for  the  coefficients  which 
appear  in  the  definition  of  Sn,i  in  terms  of  the  Fourier  coefficients  of  /.  Several  results 
concerning  the  behavior  and  convergence  of  the  approximations  are  presented  in 

section  4.  In  sections  5-6,  the  results  of  sections  2-4  are  generalized  to  values  of  Af  >  2. 
In  section  7  the  approximations  are  applied  to  a  simple  initial,  boundary-value  problem 
for  the  heat  equation.  Some  observations  and  insights  about  this  class  of  approximations 
are  discussed  in  section  8. 
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2.  Odd  functions  of  x.  Since  /  can  be  expressed  as  the  sum  of  an  even  function 
of  X  and  an  odd  function  of  x,  it  is  sufficient  to  consider  separately  the  cases  when  / 
is  either  even  or  odd.  Thus,  in  this  section,  we  consider  the  case  when  /  is  an  odd 
function  of  x,  while  the  case  of  an  even  function  will  be  discussed  in  the  next  section. 
For  this  case,  the  Fourier  series  of  /  can  be  expressed  as  (1)  with  each  a„  =  0,  n  >  0, 
and  we  define  a  class  of  Fourier-Pade  approximations  SsMi^)  of  form 


(7) 


c  _  EiLi  Bn  8in(nx) 

1  +  Cm  cos(mx)  ■ 


In  this  section  we  examine  in  detail  the  case  Af  =  1.  The  case  M  >2  will  be  discussed 
in  section  5.  Thus  we  consider  first  approximations  of  the  form 


(8) 


Sn,i{x) 


ESLi^nSin(nx) 
1  +  Cl  cos(x) 


The  constants  are  determined  by  conditions  (6),  which  for  this 

case  become 


(9) 


where 


N 


Y^dkjBj  =  bk, 

3=1 


fc  =  l, 


N+1, 


dkj  =  (1/^r) 


sin(jx)8in(fcx) 
1  +  Cl  cos(x) 


dx  =  dj,k. 


To  solve  (9),  we  first  use  the  identity 


sin(jx)  sin(A:x)  =  (l/2)(cos((A:  —  j)x)  —  cos((A:  +  i)x) 
to  express  the  coefficients  {dfcj}  as 

(10)  dkj  =  /|fc-j|,i  -  /fc+i.i,  where  /„,i  =  (I/tt)  ^  ^  0- 

Assuming  that  0  <  C^  <  1,  we  find  (see,  e.g.,  [8],  p.  113) 

(11)  /„,i  =  ,  for  n  >  0,  where  p  =  - -  or  Ci  =  -  ^  . 

I  —  p*  Oi  1  +  p* 

Using  this  result  in  (10)  we  can  write 

dk^  =  /"■'(I  +  p^)(l  +P^  +  for  k  >  j. 

Using  these  expressions,  (9)  can  be  expressed  in  matrix-vector  fotm  as 


( 

(  ^ 

(12) 

% 

Ba 

1 

b2 

k  Bn  > 

~  1+P^ 

<  bff+i  > 

3 


where  the  {N  +  1)  x  N  matrix  R  is  defined  by 


/I  p 

p  1  +  p(l  +  p2) 

p(l  +  p^)  1  +  p^  +  p^ 


p^-» 

p^-»(l+p^) 
p^-®(l  +  p2  +  p^) 


p^-l  p^-2(l+^2) 

V  P^  P^-HI  +  P^) 


l+p2^...^p2N-2 

p(l  +  p*  +  •  •  •  +  / 


To  solve  (12),  we  first  note  that  if  we  subtract  p  times  the  next  to  last  equation  from 
the  last  equation  we  obtain  the  relation  0  =  (6w+i  —  />^Jv)/(l  +  p^),  from  which  we  find 
(assuming  ^  0) 


—  r  —  —  ^^Nbs+i 

^  bs  '  *  1  +  p*  bff^  + 

In  a  similar  manner,  subtracting  p  times  the  (k  —  1)‘*  equation  from  the  equation 
we  find 


bk  -  pbk-i 

1  +  p2 


N 


-  E 

i=fc+i 


1  <  Jb  <  Af. 


Using  this  recursive  definition  of  the  {Bk},  we  can  use  induction  to  show  that  the  {Bk} 
are  given  explicitly  by 

(14)  Bi.=K-  =l>k+  ^(^-1  +  6k+i)  1  <  t  <  at, 

1  +P^  2 

(We  note  that  equations  (14)  also  follow  from  equating  (8)  to  sin(jx),  multiply¬ 

ing  by  1  -f  Cl  cos(x),  and  equating  coefficients  of  sin  (nx)  on  each  side  of  the  resulting 
expression,  for  1  <  n  <  N.)  Thus,  the  coefficients  {Bjt}  and  Ci,  which  appear  in 
the  definition  (8)  of  Ss,i,  are  completely  defined  in  terms  of  the  coefficients  {bk}  by 
equations  (13)  and  (14). 

Example  1:  As  an  example  of  these  results,  we  let  /(x)  =  x/2,  for  — tt  <  x  <  x, 
and  /(x  2x)  =  /(x).  Then  6*  =  (— l)*+*/fc,  for  k  =  1,2,...  .  Using  the  formulas 

above,  we  find  p  =  —N/{N  +  1)  and  hence 

2N{N-kl)  3N{N  +  l)  +  2 

'  2N(N  +  1)  +  1'  ^  AN{N-\-l)  +  2' 


Bk  =  (-1)* 


2N{N  -1-1)  jb 
2Ar(W -I- 1) 1  Jb2  -  1 


,  A:  =  2, 3, . . . ,  AT. 


Explicit  values  for  Ci  and  Bi, ...,  Bn  for  a  few  small  values  of  N  are  shown  in  Table  1. 
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TABLE  I  {Example  1:  M  =  1) 


N 

Cl 

Bx 

B3 

B3 

Ba 

B3 

Be 

1 

4/5 

4/5 

2 

12/13 

10/13 

3/26 

3 

24/25 

19/25 

7/50 

-2/75 

4 

40/41 

31/41 

37/246 

-4/123 

5/492 

5 

60/61 

46/61 

19/122 

-13/366 

3/244 

-3/610 

6 

84/85 

64/85 

21  im 

-19/510 

23/1700 

-1/170 

7/2550 

To  illustrate  these  results,  in  Figure  1  we  have  plotted  the  function  /(x)  (solid  line), 
along  with  the  approximations  54,o(x)  (dashed  line)  and  53,i(x)  (dotted  line).  As  the 
figure  illustrates,  except  possibly  in  a  neighborhood  of  the  discontinuities  of  /,  53,i  is 
a  better  approximation  to  /  than  is  54,0.  In  particular,  the  oscillations  present  in  54,0 
away  from  the  singularities  of  /  have  been  virtually  eliminated  in  83,1.  However,  near 
the  points  of  discontinuity  of  /,  the  Gibbs  phenomenon  is  still  present  in  534,  although 
its  effect  appears  to  be  somewhat  less  than  the  corresponding  phenomena  present  in 
84,0.  We  shall  discuss  this  example  further  in  sections  4  and  8. 

Before  continuing,  we  briefly  note  two  special  cases.  First,  if  bs+i  =  0,  then  p  =  0 
and  hence  Ci  =  0.  Thus  Bk  =  bk,  for  1  <  fc  <  AT,  and  8s, i  reduces  to  the  Fourier  partial 
sum  8s, 0,  for  this  case.  Second,  suppose  6*  =  0  when  k  is  even.  In  this  case,  we  modify 
the  form  of  8s,i  and  define 


(15) 


Ss,i{x) 


T,n=l  B2n-l  8in((2n  -  l)x) 
1  +  C2  C08(2x) 


Then  equations  (6)  for  this  case  have  the  same  form  as  (9)  with  bk  replaced  by  63;^-! 
and  dkj  =  dj^k  =  I\k-j\  —  h+j-i^  from  which  we  find 


fo,  k  >  j. 

1 

Using  these  expressions  and  the  same  type  of  arguments  as  presented  above,  we  find 
for  this  case  (assuming  b^s-x  0) 


(16)  C-i  —  —  j  ^ ^ ,  p  =  +  ~^{^k-3  +  frjib+l))  1  <k  <  N. 

(Recall  that  we  have  defined  6_i  =  —61.) 
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3.  Even  functions  of  x.  When  /  is  an  even  function  of  x,  its  Fourier  series  has  the 
form  (1)  with  each  6n  =  0.  Then  we  define  the  family  of  Fourier- Fade  approximations 


(17) 


^  ,  A)/2-l-En=l^nCOs(nx) 


In  this  section  we  consider  the  case  M  =  1.  The  case  Af  >  2  will  be  discussed  in  section 
5.  Thus  we  consider  approximations  of  the  form 


(18) 


^  /  X  >lo/2  +  EilLl  cos(nx) 

= - l'+c7^s{7) 


The  requirement  that  the  first  -|-  2  Fourier  coefficients  of  Sn,i  agree  with  those  of  / 
leads  to  a  system  of  equations  for  the  TV  -f-  2  unknowns  {Ci,Ao, >4^},  which  can  be 
solved  using  the  techniques  of  section  2.  In  particular,  we  find  (assuming  as  ^  0) 


(19)  Cl  =  -  where  p  =  Ak  =  ak  +  ^(<**-1  +  a*+i),  0<k<N. 

l-r  p*  oln  ^ 

In  the  special  case  when  nv+i  =  0,  we  see  that  p  =  0  and  hence  C\  —  0.  Then  each 
Ak  =  flfc  and  Ss.i{x)  =  Ssjai^)  for  this  case.  Also,  in  the  special  case  that  ajk  =  0  when 
k  is  odd,  we  redefine 


(20) 


C  (  X  +  EiLi  cos(2nx) 

*  l-bC2COS(2x) 


and  find  (assuming  a2s  ^  0) 

(21)C'3  =  -  where  p  =  y42*  =  02*  +  §(a2*-2  +  02^+2),  0<k<N. 

P*  «27V  2 

Example  2:  As  an  application  of  these  results,  we  let  /(x)  =  1/^1  -I-  aco8(x), 
where  <  1.  Here  /  is  analytic  for  — ir  <  x  <  x,  but  develops  a  sharp,  narrow  peak 
near  x  =  0  as  a  — ►  —1.  In  Table  II  we  have  recorded  the  Fourier  coefficients  {on}  of  /, 
as  well  as  the  coefficients  {C\,  Ao, ...,  for  a  few  small  values  of  N  when  a  =  —0.95. 


Table  II  {Example  2:  o  =  —0.95,  Af  =  1) 


1  ^ 

as 

Cl 

Ao 

A, 

A2 

A3 

A4 

2.94734 

-0.68380 

2.15059 

1 

1.16518 

-0.85289 

1.95357 

-0.37012 

2 

0.65289 

-0.89142 

1.90867 

-0.43948 

-0.04495 

3 

0.40050 

-0.90812 

1.88922 

-0.46953 

-0.05802 

-0.01235 

4 

0.25636 

-0.91741 

1.87839 

-0.48626 

-0.06529 

-0.01657 

-0.00452 
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To  illustrate  these  results,  in  Figure  2  we  have  plotted  the  function  f{x)  (solid 
line),  along  with  the  approximations  S4ft{x)  (dashed  line)  and  53,i(x)  (dotted  line).  As 
the  figure  illustrates,  53,1  is  a  noticeably  better  approximation  to  /  than  is  84,0.  In 
particular,  the  oscillations  present  in  84,0  have  been  virtually  eliminated  in  83,1.  We 
shall  discuss  this  example  further  in  sections  4  and  8. 

4.  Analysis  of  the  Case  M  =  1.  We  now  prove  certain  results  which  express 
the  asymptotic  behavior  of  the  coefficients  Ci  and  {Afc}  or  as  N  —*  00,  and  which 
also  allow  us  to  make  some  statements  concerning  the  convergence  of  the  family  of 
approximations  {Stv.i},  as  N  00. 

It  is  well  known  (see  [2],  for  example)  that  if  /  has  q  continuous  derivatives  on 
—w<x<ir  and  the  derivative  of  /  of  order  9  +  1  is  piece- wise  continuous,  then  the 
Fourier  coefficients  of  /  are  0(1/Jfc®+*)  as  k  —*  00.  For  such  a  function,  we  now  show 
that  the  coefficients  in  8s,i  decay  as  Jb  — »  00  at  a  faster  rate  than  the  Fourier  coefficients 
on  which  they  are  based. 

Theorem  1.  Let  f  be  an  odd,  2ir -periodic,  piece-wise  smooth  function  and  let  its 
Fourier  coefficients  {6*}  satisfy  the  condition  that  6*  =  0(1/ k^),  as  k  00,  where  p  is 
a  real,  positive  number.  Then  (Ci|  =  1-1-  0(1/Ar*)  and  Bk  =  0{l/k^^)  as  k,N  —*  00. 
More  precisely,  let  b,  $,  and  7  be  constants  (independent  ofk).  Then,  ifbk  =  (6/ <?**)(!  -f 
Pjk  + ')lk'^  0(1/^)),  as  k  00,  it  follows  that 

as  k,  N  —*  00.  If  f  is  an  even,  2ir-periodic,  piece-wise  smooth  function,  then  these  results 
hold  with  bk  replaced  by  ak  and  Bk  replaced  by  Ak. 

Proof:  We  shall  outline  the  proof  only  for  the  case  when  /  is  an  odd,  27r-periodic, 
piece-wise  smooth  function,  since  the  proof  for  the  case  when  /  is  an  even  function 
follows  the  same  line  of  reasoning.  Using  the  assumed  form  of  the  asymptotic  behavior 
of  the  coefficients  {bk},  we  note  first  that  we  can  write 

6w+i  =  W(N  +  lY)(l  +  ff/(N  +  l)  +  y/(N-Hf  +  0{l/{N  +  lf)) 

(22)  =  (6/JV>’)(l+e)->'(l +/?€/(! +  e)  +  7£V(l+')'  +  0(l/JV’)) 

=  (6/JV)(l  +  (^  -  p)e  +  (7  -  /?(p+  1)  +  p(p  +  l)/2)e‘  +  0(1/JV»)), 

where  we  have  defined  e  =  l/N.  Then,  using  this  expression  and  the  definition  of  p  in 
(13),  we  can  write 

(23)  p=‘^  =  l_£  +  E±Z^  +  0(l/jV^),„jV^oo, 
and  hence 

as  stated  in  the  Theorem. 
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In  a  similar  manner,  for  large  values  of  k  we  can  use  the  assumed  form  of  the 
coefficients  {6^}  to  write 

bk±i  =  mk±iy)ii+mk±i)+i/{k±iy+o{i/k^)) 

=  {b/l^)i\  ±  €)-Pil  +  0€/il  i  e)  +  7eV(l  ±  +  Oil/l^)) 

=  ib/l^)(l  +  (^  :f  p)e  +  (7  q:  +  1)  +  p(p  +  l)/2y  +  0{\/k^)), 

where  we  have  now  defined  e  =  l/k.  Using  these  expressions,  along  with  an  analogous 
expression  for  bk  and  the  asymptotic  expression  above  for  Ci,  in  the  definition  of  Bk  in 
(14),  we  find  for  large  values  of  k  and  N  (with  k  <  N) 


Bk  =  bk  +  (C\l2){bk~x  +  6*+i) 

=  {b/kP){l  +  ^t  +  ^t^  +  Oi^) 

+(l/2)(-l  +  p"/(2N^)  +  0(1/A^3))(2  +  2$t  +  (27  +  p  +  p^)e‘^  +  0(c3))} 

=  -(6p/2)(l/fc»-+2)  (1  +  p(i  _  ^ 

which  completes  our  proof. 

Thus,  Theorem  1  es^blishes  that  the  coefficients  {^n}  in  the  approximations  Ss,\ 
decay  to  zero  more  rapidly  (by  a  factor  of  k~'^)  than  the  coefficients  in  the  partial  sum 
Sn,q  on  which  they  are  based. 

Before  continuing,  we  remark  that  results  very  similar  to  those  stated  in  Theorem 
1  hold  if  the  coefficients  {5*}  of  /  have  the  same  asymptotic  form  as  that  stated  in  the 
theorem,  except  for  a  multiplicative  factor  of  (—1)*,  i.e.,  if  bk  =  (6/^*’)(— 1)*(1  +  /?/A  + 
“1 ! 0{\ ! l^)),  dA  k  oo.  In  this  case,  it  is  easy  to  show  that  the  Fourier  coefficients 
(bk)  of  /  =  f{x  +  tt)  are  related  to  those  of  /  by  5*  =  (— 1)*5*.  It  then  follows  that 
the  coefficients  {6*}  satisfy  the  conditions  of  the  theorem  and  hence  the  corresponding 
coefficients  C\  and  {Bk}  have  the  asymptotic  form  indicated  in  the  theorem.  In  partic¬ 
ular,  it  follows  from  equations  (13)  and  (14)  that  Ci  =  —Ci  and  Bk  =  (— 1)*5*.  Thus, 
the  asymptotic  behavior  of  Ci  is  the  negative  of  that  indicated  for  Cj  in  the  theorem, 
while  the  asymptotic  form  of  Bk  is  (—1)*  times  the  form  indicated  for  Bk- 

To  illustrate  these  results,  in  Figure  3  we  have  used  the  coefficients  {Bfc}  from 
Example  1  and  have  plotted  k^Bk\  as  a  function  of  1/A:,  for  3  <  iV  <  20.  The  figure 
clearly  illustrates  the  decay  of  the  coefficients  {Bk},  and  also  indicates  an 

interesting  asymptotic  behavior  of  Bjv_„  for  a  fixed  value  of  g,  as  JV  -♦  oo. 

If  the  function  /  is  analytic  for  — jt  <  x  <  tt,  then  (see,  e.g.,  [2])  the  Fourier 
coefficients  of  /  decay  exponentially  as  k  oo.  That  is,  there  exists  a  constant  0,  with 
0  <  1^1  <  1,  such  that  are  0(^*)  as  k  oo.  The  next  theorem  shows  that, 

again  for  this  case,  the  coefficients  {Bk}  decay  to  zero  more  rapidly  than  the  Fourier 
coefficients  on  which  they  are  based. 

Theorem  2.  Let  f  be  an  odd,  analytic,  2n -periodic  function  and  let  its  Fourier 
coefficients  {bk}  satisfy  the  condition  that  bk  =  0{\9\’‘/kP),  ask  oo,  where  0  <  |^(  <  1, 
and  p  is  a  real  number.  Then  Ci  =  20 /{I  +  9^)  -h  0{lfN)  and  Bk  =  0{\9\^  j as 
k,N  -*  oo.  More  precisely,  let  b,  fi,  and  7  be  constants  (independent  of  k).  Then,  if 
bk  =  ib0'‘lkP){l  +  j3/k  -I-  7/P  -f  0(1/P)),  as  k  —*  cx>,  it  follows  that 


20  f  p  (p-2/3)(l-g^)-^pHl-6P-^g^) 

N\l-\-0^)  2(1 -I- 
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Bk  = 


“  1  \\r^Wp+{20  +  p)(l  +  k/Nm  -  k/N) 


-2(IW  (<* + + ^>)  0  -  ‘/Af)  -  rfji] } + 0(;^). 

as  k,N  —*  oo.  If  f  is  an  even,  analytic,  2v-periodic  function,  then  these  results  hold 
with  bk  replaced  by  ak  and  Bk  replaced  by  Ak- 

Proof:  As  in  the  proof  of  Theorem  1,  we  shall  outline  the  proof  only  for  the  case 
when  /  is  an  odd  function,  since  the  proof  for  the  case  when  /  is  an  even  function 
follows  the  same  line  of  reasoning.  Using  the  assumed  form  of  the  asymptotic  behavior 
of  the  coefficients  we  note  first  that  bii/^i  h2is  the  asymptotic  form  indicated  on 
the  right  side  of  equation  (22),  multiplied  by  6^'*’^.  Then,  using  this  expression  in  the 
definition  of  p  in  (13),  we  find  that  p  has  the  asymptotic  form  indicated  by  the  right 
side  of  (23)  multiplied  by  6.  Using  these  expressions  in  the  definition  of  Ci  in  (13)  we 
find 

p _ 2p _  l-p/N  +  (p  +  p‘-2m2N^)  +  0{im 

'  l+P^  l+em-p/N  +  {p  +  ^ -20)/i,2N^)  +  O(l/N^)f 

_ 2»_  (,_p_  .  (p-2m-e')+p‘(i-ie^  +  e')  .  f,,,  ,^3,1 

l+«2\  JVV1+9V'''  2(1+«S)W 

as  N  —*  oo. 

In  a  similar  manner,  using  the  assumed  asymptotic  form  of  the  coefficients  {6^}  we 
can  write 

bka  =  +  (^  TP)e  +  (7  T  ;8(P  +  1)  +  P{P  +  l)/2)£=  +  0(1/A:’)), 

as  k  —b  OO,  where  c  =  l/k.  Thus,  for  large  values  of  both  k  and  N  we  have 

Bk  =  bk  +  {Cif2){bk-i  +  6jfe+i) 


$^b(  ^  ,  if  p  Ii-e^\ 

~  kp  “rT^P“]v(rT^)'^ 


(p-2^)(l-g^)-fp"(l-6g"+g") 

2(1+02)2^2 


•  [l  +  (^  +  p)c  +  (7  +  ^{p  +  1)  +  p{p  +  l)/2)c* 

+  02  (l  +  (^  -  p)e  +  (7  -  0{p  +  1)  +  p{p  +  l)/2)c2)] }  +  O(0Vfc^") 


=  {'’IT?**  -  -  I  +  (2/3  +  p)(l  +  (:/JV))(l  -  k/N) 


-2(rw  -«'’+«*))(!-  kIN)  -  ^] } + O(^), 

ds  k  —*  oo.  This  completes  our  proof. 

We  note  that  the  rate  at  which  the  coefficients  {Bjt}  decay  to  zero  changes  slightly 
as  A:  — »  W.  In  particular,  for  k  <  Bk  ~  0($^(1  —  0^)(1  —  klN)/k'^^y  but  Bs  = 
O{0^ fN^^),  as  k,N  —*  oo.  Thus,  while  each  Bk  decays  to  zero  f^lster  than  the  corre¬ 
sponding  Fourier  coefficient  6*,  we  see  that  decays  somewhat  faster  than  B*,,  when 
k2  <  ki  <  N. 

To  illustrate  these  results,  in  Figure  4  we  have  used  the  coefficients  {Afc}  from 
Example  2  and  have  plotted  A:^lAfc|/|aifc|  as  a  function  of  1/ A:,  for  3  <  N  <  20.  The 
figure  clearly  illustrates  that  the  coefficients  {y4fc}  decay  at  a  faster  rate  (asymptotically 
by  a  factor  of  A:"^)  than  the  coefficients  {a*},  and  also  indicates,  as  in  Example  1,  an 
interesting  asymptotic  behavior  of  A^-q,  for  a  fixed  value  of  q,  as  N  —*  oo. 

Finally,  we  use  the  results  of  the  previous  two  theorems  to  show  the  manner  in 
which  the  approximations  {5Ar,i(x)}  converge  to  the  original  function  /(x). 

Theorem  3.  Let  f  be  an  odd,  2ir-periodic,  piece-wise  smooth  function  and  let 
its  Fourier  coefficients  {6*}  satisfy  the  conditions  of  Theorem  1.  Then  the  sequence  of 
approximations  {5a^,i(x)}  converges  to  (/(x"*")  -1-  /(x~))/2,  as  N  —*  oo,  for  all  — t  < 
X  <  IT.  Moreover,  the  sequence  {5a^,i(x)}  converges  like  a  Fourier  series  whose  terms 
are  0(l/iVP+^).  These  results  also  hold  if  f  is  an  even,  2‘K-periodic,  piece-wise  smooth 
function,  if  the  Fourier  coefficients  {a*}  of  f  satisfy  the  conditions  stated  above  with  bk 
replaced  by  Ok. 

If  f  is  an  odd,  analytic,  2v -periodic  function  and  its  Fourier  coefficients  satisfy  the 
conditions  of  Theorem  2,  then  the  sequence  of  approximations  {■S'a^.iCx)}  converges  to 
/(x),  as  N  —k  oo,  for  all  — tt  <  x  <  tt.  Moreover,  the  sequence  {5;v,i(®)}  converges  like 
a  Fourier  series  whose  terms  are  O{0^  These  results  also  hold  if  f  is  an  even, 

analytic,  2ir-periodic  function,  if  the  Fourier  coefficients  {cfc}  of  f  satisfy  the  conditions 
stated  above  with  bk  replaced  by  ak. 

Proof.  As  in  the  proofs  of  the  previous  theorems,  we  shall  outline  only  the  proof 
for  the  case  when  /  is  an  odd,  27r-periodic,  piece-wise  smooth  function. 

To  begin,  we  note  that,  since  /  is  an  odd  function,  we  can  restrict  our  attention  to 
values  of  x  in  the  open  interval  0  <  x  <  ^  (since  5a^,i(0)  =  (/(O'*")  +  /(0“))/2  =  0  = 
5'Ar,i(’r)  =  (/(tt'*')  +  f{ir~))l2,  and  hence  convergence  is  assured  at  x  =  0  and  x  =  tt  in 
a  trivial  way).  Consequently,  we  let  x  be  any  fixed  number  in  the  interval  0  <  x  <  tt 
and  let  e  >  0  be  any  fixed  positive  number.  We  then  define 

(24)  =  {fix*)  +  f{x-))/2  -  S«..(a:)  =  Bg'Cx)  +  E*^\x), 

(25)  EH'ix)  =  (fix*)  +  fix-))/2  -  S„,„(x),  f:g'(x)  =  SKfl(x)  -  S«,.(i). 

By  the  hypotheses  of  the  theorem  concerning  /,  the  Fourier  partial  sums  Sn,o{x) 
converge  to  (/(x+)  -I-  /(x“))/2,  as  N  -*  oo.  Hence,  there  exists  a  positive  integer 
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M  =  Ni{x,t)  such  that 


=  |(/(x+)  +  f(x-))/2  -  5jv.o(x)|  <  for  all  N  >  Ni. 


Also,  using  the  definitions  of  5w,o  5^,1,  we  can  write 


eS'M  =  S«.o(i)  -  =  E  ».8m(nx)  - 

1  +  Cl  cos(x) 

/27\  _  Cl  bs  sm({N  +  l)ar)  -  bj^+i  sinjNx) 

^  ^  2  l  +  CiCos(x) 

Using  the  assumed  asymptotic  form  of  the  coefficients  {6n}»  there  exists  a  positive 
integer  IVj  =  N2ix,  e)  >  (4 16|  /(€(1  —  cos(ar)))^^’’  such  that  |6nj  <  2  \b\  /N^  and  Ci  >  —1 
for  all  JV  >  A/j.  Then,  using  equation  (27),  we  can  write 

|4“'(x)|  =  |5«4x)  -  S«,.(x)|  < 

Using  the  bounds  (26)  and  (28),  from  equation  (25)  we  find 

(29)  |(/(x*)  +  /(i-))/2  -  5«.,(x)|  <  |£:<;>(x)|  +  |£:jj>(x)|  <  |  +  |  =  t, 

for  all  TV  >  A/3  =  Max{Ni,  N2).  Thus,  the  sequence  of  approximations  {5^,i}  converges 
to  (/(x"*")  +  /(x~))/2  as  A/  — »  00. 

To  demonstrate  the  manner  in  which  these  approximations  converge,  we  write  Sn,i 


4  V 

(30)  Sn,i  =  5^  r*,  where  n  =  5i,i  and  r*  =  5fc,i(x)  -  5fe_i,i(x),  2<k<N, 

k=l 

and  then  examine  the  quantities  r*  for  large  values  of  k.  Using  the  definition  (8),  we 
can  write 


rifc(x)  = 


E*=i  Bi>‘)  sin(nx)  8in(nx) 

l  +  C{*^cos(ar)  1  +  cos(x) 


where  we  have  placed  a  superscript  on  the  coefficients  Bn  and  Ci,  to  remind  us  of  the 
dependence  of  these  quantities  on  the  parameter  k.  Combining  the  terms  on  the  right 
side  of  equation  (31)  and  using  various  trigonometric  identities,  we  can  write 


rfe(x)  = 


Num 
Den  ’ 


where 


Num  =  ”»((*- »)^)ci‘) 


aiB(kx) 


{^(i-1)  Mk)  'j 

c!‘-‘>|6t  +  ^(»,.,+kw.)|, 


sin((fc+  l)x) 

2 

and 

(33)  Den  =  (1  +  C{*>  cos(x))(l  +  co8(x)). 

Using  the  asymptotic  form  of  the  coefficients  and  C|*^  from  Theorem  1,  we  find 


Num  = 


bp  1 
'~2k^ 


(1  —  cos(x))  sin(A:a;)  + 


(34)  Den  =  (1  —  cos(x))^  +  ^(1  —  cos(x))  cos(x)  +  0(l/fc^),  as  k  -*  oo. 

k* 

Using  the  estimates  (34)  in  (32),  we  find 

This  completes  the  proof  of  our  theorem. 

To  illustrate  some  of  these  results,  in  Figure  5  we  have  plotted  loglD;v,i|/log|^l 
for  Example  1,  as  a  function  of  1/  log  \N\,  {ot3<  N  <  50,  with  x  =  7r/2.  We  note  that, 
if  En,i  ~  aN~^,  as  AT  -♦  oo,  then  log  |DAf,il/log  \N\  ~  -p  +  log  la]/ log  |iV|,  as  N  -*  oo. 
Thus,  the  plot  shotild  be  approximately  linear  for  small  values  of  l/log|7V|,  with  an 
intercept  of  —p.  x  he  value  of  p  =  3,  as  predicted  by  Theorem  3,  is  clearly  suggested  by 
the  figure.  In  this  figure  we  have  also  plotted  the  corresponding  error  for  the  Fourier 
partial  sums,  corresponding  to  A/  =  0.  For  M  =  0,  the  figure  suggests  a  value  of  p  =  1. 
Thus,  the  improvement  in  the  rate  of  convergence  suggested  by  the  figure  is  consistent 
with  the  increase  predicted  by  Theorem  3,  i.e.,  an  increase  in  p  of  2. 

5.  Case  M  >2.  We  now  generalize  the  results  of  the  previous  sections.  In  par¬ 
ticular,  when  /  is  an  odd  function  of  x,  we  consider  the  family  of  approximations 


Sn,m{x) 


Enr:!  sin(nx) 

1  +i;m=lC^mCOs(mx)’ 
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(36) 


where  M  is  any  positive  integer.  The  constants  {Bn^Cm}  are  determined  by  the  con¬ 
ditions  (6),  which  for  this  case  become 


(37) 

where 


s 


^  =  1>  *1 N  +  M, 

1=1 


(38)  =  (IM  - 


'  sin(jx)sin(Ara:) 


+  Em=l  C08(mx) 


dx  =  dyk  - 


(39) 


r  /  X  cos(nx)  , 

In,M  =  (I/’t)  /  - — rnj — ^ - ;; - rdx,  n  >  0. 


To  evaluate  /n,M,  we  denote  the  denominator  of  the  integrand  by  P  and  express  cos(mx) 
as  a  polynomial  of  degree  m  in  cos(x).  Then  P  becomes  a  polynomial  of  degree  M  in 
cos(x)  and,  aissuming  Cm  ^  0,  we  can  write 


(40) 

and  hence 

(41) 


M  M 

P  =  l+'^Cm  cos(mx)  =  ^~^Cm  n(«>s(x)  -h  X,), 

msl  «ssl 


M 


M 


1  +  Em=i  Cm  cos(mx)  2^-^Cm  cos(x)  -I-  Zi 


kslfk^i 


Here  {— are  the  roots  of  the  polynomial  P,  when  P  is  regarded  as  a  polynomial  in  the 
variable  z  =  cos(x),  and  we  have  assumed  that  the  {z,}  are  all  distinct.  Substituting 
(41)  into  (39)  and  using  (11),  we  can  write 


(42) 


9  ^ 

*  «  Pi 


T  —  ^  'ST'  a  t'i 

dn,M  ^  /  .  Pi  1  2  > 

Cm  S  P  i 


M 

n 


PiPk 


iPk  -  Pi)i^  ~  PiPk) 
Using  this  result  in  (38)  we  find 


,  Pi  =  \Jzf  -  1  -  Zi. 


(43) 


*>j- 


1=0 


To  solve  equations  (37),  we  first  define  the  quantities 

5o  =  1,  Sm  =  (-1)"*  2^  Pi,  •  •  ■  Ptm  ,  1  <  m  <  M. 

l<il  <i2<...<im<Ar 
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(44) 


(We  note  that  is  just  the  coefficient  of  in  the  polynomial  ni^i(p  ~  Pi)-)  We 
then  denote  the  equation  in  (37)  by  eq{k)  and  form  the  linear  combinations 

M 

(45)  lc{k)  =  s^eqik  -  m),  k  =  N  M,N  +  M  -  I, 1. 

m=0 

Here  we  define  69(0)  to  be  the  trivial  equation  0  =  0,  and  eq{n)  =  —eq{—n),  for  n  <  0. 
For  k  =  N  +  M,N-\-M  —  1, IV  +  1,  we  find  that  the  left  side  of  each  lc{k)  vanishes. 
This  leads  to  the  relations 
M 

(46)  ^m^TV+M+l-p-m  =  — hN+M+l-p>  P  =  1, 

m=l 

which  are  a  system  of  M  linear  equations  for  the  M  quantities  We  shall 

assume  that  the  determinant  of  the  matrix  of  coefficients  in  (46)  is  nonzero  so  that 
these  equations  can  be  solved  uniquely  for  Once  these  quantities  have  been 

determined,  we  can  use  equations  (44)  and  (42)  to  express  each  s,n  in  terms  of  the 
quantities  z,,  and  then  use  equation  (41)  to  express  each  of  the  coefficients  Cm  in  terms 
of  the  Zj.  We  write  out  these  equations  for  a  few  small  values  of  M  and  then  use  induction 
to  show  that  the  coefficients  {Cm}  can  be  expressed  explicitly  in  terms  of  the  known 
quantities  {s,-}  as 

(47)  =  4  E  l<m<M,  D  =  Y2^]- 

^  k=0  3=0 

To  determine  the  coefficients  {H*},  for  k  =  N^N  —  we  solve  lc{k)  for  Bk 

and  then  use  induction  to  show  that  the  coefficients  {Bk)  are  given  explicitly  by 

(48)  Bk  =  hk  +  \^  Cm{bk+m  +  bk-m),  I  <  k  <  N. 

(Equations  (48)  also  follow  from  equating  (36)  to  bj  sin(jx),  multiplying  by  1  + 

53m=i  Cm  cos(mx),  and  equating  coefficients  of  sin  (nx)  on  each  side  of  the  resiilting 
expression,  for  1  <  n  <  N.)  Thus,  the  coefficients  {Cm}  and  {Bn}  are  completely 
determined  by  equations  (48),  (47),  and  (46). 

Example  3:  As  an  application  of  these  results,  we  consider  the  odd,  2x— periodic 
function  /  defined  by  /(x)  =  ex(x  —  x)f{€^  +  (x  —  5)*(x  —  x  +  6)^)^/^  for  0  <  x  <  x, 
where  c  and  6  are  real,  positive  parameters.  For  this  function,  we  find  ft*  =  0  when  k  is 
even  and  hence  Sn,o  =  sin((2fc  —  l)x).  Then,  for  Af  =  2,  we  define 

E  B2k-i  8in((2*:  -  l)x) 
c  _  fc=i _ 

1  +  Cl  cos(2x)  +  C2  co8(4x)  ’ 

where  the  coefficients  [B^k-i)  and  {Cm}  are  defined  by  equations  (46)-(48),  with  bj 
replaced  by  62,-1.  Explicit  values  for  {Ci,C2, Bi, B2N-i}  are  shown  for  a  few  small 
values  of  N  in  Table  III  when  c  =  0.1  and  6  =  0.5. 
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Table  III  {M  =  2) 


N 

biN-X 

Cl 

C2 

Bi 

^3 

Bs 

B-r 

Bit 

1 

0.41814 

-.99477 

0.29067 

0.43574 

2 

0.31809 

-1.3981 

0.46568 

0.43655 

-.16521 

3 

0.09684 

-1.3559 

0.54610 

0.42554 

-.17564 

0.02412 

4 

-.11150 

-1.3462 

0.58487 

0.42077 

-.18343 

0.02989 

0.00716 

5 

-.17153 

-1.3539 

0.59107 

0.42047 

-.18703 

0.02987 

0.00817 

0.00234 

In  Figure  6,  we  have  plotted  56, o,  <^5.1  (using  equation  (15)),  and  54, j,  as  well  as  /,  for 
0  <  X  <  IT.  From  the  figure  it  is  clear  that  54,2  is  a  better  approximation  to  /  than 
either  55,1  or  56,o.  In  fact,  for  this  example,  5s, 1  and  56,o  lie  close  to  each  other  and 
are  each  a  ‘‘poor”  approximation  to  /.  We  shall  comment  further  on  this  example  in 
section  8. 

If  /  is  an  even  function  of  x,  then  results  very  similar  to  those  above  hold  with  each 
bje  replaced  by  at..  In  particular,  for  approximations  of  the  form 


we  find 


Snm(^) 


A)/2  +  EiLi>^nC08(nx) 

1  + 1^=1  Cm  cos(mx) 


1  " 

=  a*  +  -  Cm(o*+m  +  a*-m)»  0  <  A:  < 

^  m=l 


o  M—m  M 

Cm  =  53  SfcSm+fc,  1  <  m  <  Af,  D  =  '^s], 

^  *=0  j=0 

where  the  {s^}  are  now  determined  from  the  equations 
M 

(49)  ^  >  SmOiV+M+l— p— m  —  OAT+M+I— p>  P  “  1,...,  A/. 

rr~\ 

6.  Analysis  of  the  case  M  >2.  In  this  section,  we  prove  some  results  analogous 
to  those  proved  in  section  4  for  the  case  Af  =  1. 

Theorem  4.  Let  f  be  an  odd,  2‘K-periodic,  piece-wise  smooth  function  and  let  its 
Fourier  coefficients  {bk}  satisfy  the  condition  that  bk  =  0{lfk^),  ask  oo,  where  p  is  a 
real,  positive  number.  Then,  for  any  positive  integer  M  >  1,  Cm  =  Cm,M-\-0{llN^),for 
1  <  m  <  M,  where  Cm,M  ore  certain  constants  independent  of  N,  and  Bk  = 
as  k,N  oo.  More  precisely,  let  b,  and  7  be  constants  (independent  ofk).  Then  if 
bk  =  (V^’’)(l  +  z®/*;  +  +  0(1  IP)),  as  k-¥  00, 

„  (-irU-m)  f.  m^(p  +  M-I)(p  +  2(M-I)) 

”  (m-i)  ^  2{2M-1)N^ 
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6p(p  +  1)  •  •  •  (p  +  A/  -  1)(-1) 


kP+2M 


n(p+2A<-y-o(*W 


+0(1/*”"'+*), 

as  k,N  —*  oo.  Here  is  the  usual  binomial  coefficient.  If  f  is  an  even, 

2‘K-periodic,  piece-wise  smooth  function,  then  these  results  hold  vnth  bk  replaced  by 
and  Bk  replaced  by  Ak- 

Proof:  As  in  the  proof  of  Theorems  1-3,  we  shall  only  outline  the  proof  for  the  case 
when  /  is  an  odd,  2T-periodic,  piece-wise  smooth  function. 

To  see  that  the  coefficients  {Cm}  have  the  stated  form  as  AT  — »  oo,  we  first  use  the 
assumed  form  of  the  asymptotic  form  of  the  coefficients  {bk}  to  write 

,  ,  =  mN+jy)(i+mN+i)+inN+jf+o[i/tp)) 

'  '  =  6€'’(1 +(yS- ip)«  +  (7-i3(p+ l)+j’p(p+ l)/2)t'  +  0(e^)), 

where  we  have  defined  e  =  1/N.  Next,  we  examine  the  behavior  of  the  quantities  {sm} 
as  N  —*  oo.  Using  the  /tn ear  equations  (46),  we  solve  for  {sm}  explicitly  in  terms  of  the 
coefficients  {A*}  for  a  few  small  values  of  M  (see  the  Appendix).  We  then  insert  the 
expansions  (50)  into  these  explicit  expressions  for  the  {sm}  and  expand  the  resulting 
expressions  for  small  values  of  e.  In  this  way  we  obtain  an  expansion  for  each  Sm  which 
is  valid  as  c  — ♦  0,  i.e.  as  AT  -+  oo.  We  then  use  induction  to  show  that  these  expressions, 
for  any  positive  integer  M,  can  be  written  as 

for  m  =  1,2,  ...,M,  as  N  —*  oo.  Inserting  equations  (51)  into  (47)  and  expanding  the 
resulting  expressions  for  large  values  of  N,  we  find  that  each  coefficient  Cm  has  the 
asymptotic  form  indicated  in  the  statement  of  vhe  theorem. 

In  a  similar  manner,  using  the  expressions  (50)  with  N  replaced  by  k,  along  with 
the  asymptotic  form  of  the  {Cm}t  in  equations  (48),  we  .'^r^d  that  the  coefficients  {Bk} 
have  the  asymptotic  form  indicated  in  the  statement  of  the  theorem,  as  both  k,N  oo. 
This  completes  the  proof  of  our  theorem. 

To  illustrate  these  results,  in  Figure  7  we  have  plotted  log|£fc|/logJb  vs.  l/logib 
for  Example  1  with  M  =  1,2,3,  and  4,  for  N  =  20.  The  figure  clearly  illustrates  how 
much  more  rapidly  the  coefficients  Bk  decay  to  zero  (especially  with  increasing  values 
of  M)  than  the  coefficients  A*,  which  decay  only  like  l/k. 

The  results  of  Theorems  2  and  3  generalize  to  the  cese  when  Af  >  2  in  a  straightfor¬ 
ward  manner.  Using  the  results  of  Theorem  4,  the  proofs  Oi  the  following  two  theorems 
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follow  closely  the  proofs  of  Theorems  2  and  3,  respectively,  and,  hence,  we  state  them 
without  proof. 

Theorem  5.  Let  f  he  an  odd,  analytic,  2^ -periodic  function  and  let  its  Fourier 
coefficients  {ft*}  satisfy  the  condition  that  ft*  =  (ft0*/ib'’)(l  ^/k  +  0{l/k^)),  as 

k  -*  oo,  where  ft,  fi,  and  7  are  constants  (independent  of  k),  p  is  a  real  positive  number, 
and  0  <  1^1  <  1.  Then 

^M-m  02k  ( M\  (  Af  \ 

C„  =  2(-ir<r  = 

Ejlo  (?) 

and  Bk  =  0(|«|*(1  -  «’)(!  -  klN)lk’*»)  +  os  t,  JV  —  oo,  ///  is  on 

even,  analytic,  2ir-periodic  function,  then  these  results  hold  with  ft*  replaced  by  a*  and 
Bk  replaced  by  Ak. 

Theorem  6.  Let  f  be  an  odd,  2^ -periodic,  piece-wise  smooth  function  and  let 
its  Fourier  coefficients  {ft*}  satisfy  the  conditions  of  Theorem  1.  Then  the  sequence  of 
approximations  {5Af,M(®)}  converges  to  (/(x'*‘)  +  f{x~))f2,  as  N  —*  00,  for  all  —ir  < 
X  <  IT.  Moreover,  the  sequence  {5/v,m(2;)}  converges  like  a  Fourier  series  whose  terms 
are  These  results  also  hold  iff  is  an  even,  2x-periodic,  piece-wise  smooth 

function,  if  the  Fourier  coefficients  {a*}  of  f  satisfy  the  conditions  stated  above  with  bk 
replaced  by  ak- 

If  f  is  an  odd,  analytic,  2ir -periodic  function  and  its  Fourier  coefficients  satisfy  the 
conditions  of  Theorem2,  then  the  sequence  of  approximations  {Sn,m(x)}  converges  to 
/(x),  as  N  -*  00,  for  all  -ir  <  x  <ir.  Moreover,  the  sequence  {Snm(^)]  converges  like 
a  Fourier  series  whose  terms  are  O{0^  fN^^),  These  results  also  hold  if  f  is  an  even, 
analytic,  2ir -periodic  function,  if  the  Fourier  coefficients  {a*}  of  f  satisfy  the  conditions 
stated  above  with  bk  replaced  by  ak. 

To  illustrate  Theorem  6,  in  Figure  5  we  have  also  plotted  log  1^1  for 

Example  1,  with  M  =  2  and  3,  as  a  function  of  l/log|,/V|,  for  3  <  <  50,  with 

X  =  7r/2.  Here  En,m(x)  =  (/(x+)  +  /(x“))/2  — 5jv,M(ic).  Using  the  ideas  discussed  after 
Theorem  3,  the  results  presented  in  Figure  5  are  consistent  with  the  convergence  rates 
of  p  =  5  (for  M  =  2)  and  p  =  7  (M  =  3),  as  predicted  by  Theorem  6  for  the  function 
considered  in  Example  1 . 

7.  Application.  Applications  of  the  results  presented  above  to  several  classes  of 
problems  usually  solved  by  Fourier  series  alone  will  be  presented  and  discussed  else¬ 
where.  In  this  section,  we  present  an  application  to  a  simple  heat  conduction  problem 
to  illustrate  some  of  the  potential  of  the  method  to  improve  the  accuracy  of  approximate 
solutions  obtained  by  partial  sums  of  Fourier  series. 

We  consider  the  problem  of  determining  the  transient  behavior  of  the  temperature 
u{x,t)  which  satisfies  the  conditions 

(52)  Ut  =  Uxx,  for  0  <  X  <  TT  and  <  >  0, 

with 


u(0,t)  =  0  =  u(ir,t),  for  all  <  >  0,  and  u(x,0)  =  x/2,  for  0  <  x  <  tt. 
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(In  equation  (52),  the  subscripts  denote  partial  differentiation.)  Using  separation  of 
variables,  the  solution  is  found  in  a  straightforward  manner  to  be 

(53)  u(ar,<)  =  ^im  u^^’°^(x,<),  where  ^ ^ — c~"**sin(nx). 

This  solution  obviously  has  a  discontinuity  at  x  =  x  when  t  =  0. 

Using  the  formulas  of  section  2,  we  define  a  new  class  of  approximate  solutions 

by 

i“i  -“i-i  ■ 

where  the  coefficients  Bn(i)  and  6'i(t)  are  defined  by  equations  (13)  and  (14)  with  bt 
replaced  by  (— l)*'‘'^c“***/^'  Thus,  we  find 

_  _  c.  „ 

(JV+l)J  +  ArJe-<«'+»)<’  ‘  4  ’ 


\  i  2  ' 

[  ib-1  ib-l-1  jj 

In  Figure  8  we  have  plotted  (dashed  line)  and  (dotted  line),  along  with  the 
exact  solution  (solid  line)  for  t  =  0.01  and  t  =  0.05.  As  the  figure  illustrates, 
is  consistently  a  better  approximation  than  to  the  exact  solution,  and  the  quality 
of  this  approximation  improves  as  t  becomes  larger.  The  corresponding  comparison  at 
t  =  0  is  the  same  as  shown  in  Figure  1  for  0  <  x  <  x. 

8.  Conclusions  and  Discussion.  We  now  make  a  few  observations  about  the 
family  of  Fourier-Pade  approximations  discussed  above,  and  tlso  comment  on  several 
of  their  properties  that  need  further  investigation. 

First,  we  note  that,  using  only  the  "information”  contained  in  the  first  few  Fourier 
coefficients  of  a  function  /,  the  functions  appear  to  provide  new  approximations 
to  /  which  are  consistently  "better”  than  the  Fourier  partial  sums  Sjvfi  on  which  they 
are  based.  The  sense  in  which  these  approximations  are  "better”  can  be  interpreted  in 
at  least  three  different  ways.  First,  Theorems  3  and  6  show  that  these  approximations 
converge  point- wise  to  (/(x"*")  +  /(x~))/2  at  a  faster  rate  than  the  original  sequence 
^  the  parameter  JV  increases.  Secondly,  although  the  Gibbs  phenomena  is 
still  present  in  the  family  {5jv,M}4  the  amplitude  of  the  oscillations  near  a  point  of 
discontinuity  of  /  appears  to  be  mitigated,  when  compared  with  the  oscillations  present 
in  Sm,o-  Finally,  oscillations  in  5/v,o,  which  lie  outside  a  neighborhood  of  a  point  of 
discontinuity  of  /  (or  outside  a  region  of  large  curvature  of  /),  are  noticeably  damped 
in  Snm^  for  Af  >  1,  especially  as  N  increases.  These  last  two  interpretations  are 
illustrated  in  Examples  1  and  2. 

The  fact  that  explicit  expressions  were  derived  for  the  coefficients  An,  Bn  and  Cm, 
which  appear  in  the  definition  of  Ss,Mi  should  be  emphasized.  (From  a  practical  point 
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of  view,  it  should  be  noted  that  it  might  be  possible  to  develop  an  efficient,  recursive 
algorithm,  similar  to  algorithms  presented  by  Baker  [1],  to  compute  these  coefficients, 
since  the  equations  which  ultimately  determine  these  quantities  are  recursive  in  nature. 
See,  especially,  equations  (46)-(48),  as  well  as  the  Appendix.)  These  explicit  expressions 
not  only  facilitate  the  proofs  of  the  stated  theorems,  but  also  illustrate  some  of  the 
potential  use  of  these  approximations  to  construct  approximate  solutions  to  problems 
involving  differential  equations,  especially  partial  differential  equations.  In  particular, 
the  Fourier  coefficients  {a„,  6„}  may  be  functions  of  one  or  more  “other”  variables,  such 
as  f,  as  in  the  example  of  section  7.  The  application  of  the  basic  ideas  presented  here  to 
problems  involving  several  different  classes  of  partial  differential  equations  is  currently 
under  investigation  and  will  be  reported  elsewhere. 

As  far  as  the  mathematical  properties  of  the  family  are  concerned,  only  the 

most  elementary  properties  have  been  investigated  here.  Many  other  questions,  which 
are  of  both  theoretical  and  practical  importance,  should  be  addressed.  For  example,  for 
a  fixed  value  of  N+M,  which  approximation  Sn,m  is  “best”?  In  Figure  9  we  have  plotted 
S6-g,q,  for  9  =  0, 1,...,5  for  Example  1.  In  this  case,  the  figure  seems  to  suggest  that 
perhaps  83,3  is  “best”  in  some  appropriate  integral  norm,  although  the  improvement  of 
the  quality  of  the  approximation  of  53,3,  say,  over  is  not  dramatic.  For  this  example, 
the  function  has  essentially  only  one  point  of  discontinuity  (at  x  =  tt)  and  hence  Sn,i 
appears  to  provide  a  “good”  approximation.  By  contrast,  the  function  considered  in 
Example  3  has  essentially  two  points  of  large  curvature  (near  x  =  S  and  x  =  ir  —  6). 
For  this  example,  Sn,i  provided  a  poor  approximation  to  /,  while  Sn,2  yields  a  much 
improved  approximation.  More  generally,  the  “best”  Sn,m  will  undoubtedly  depend  on 
both  the  general  “shape”  and  smoothness  properties  of  the  function  /,  as  well  as  on 
the  particular  norm  used. 

From  a  practical  point  of  view,  a  related  issue  concerns  how  best  to  represent  a 
function  with  a  discontinuity  interior  to  [— tt,  tt].  For  example,  consider  the  odd,  piece- 
wise  continuous,  27r— periodic  function  /  defined  by  f{x)  =  irxf2,  for  — jr/2  <  x  <  jr/2, 
and  f{x)  —  0,  otherwise,  in  [— 7r,7r].  The  Fourier  coefficients  {6„}  of  /  are  given  by 
bn  =  (— if  n  is  odd,  and  6„  =  (x’/(2n))(— if  n  is  even.  The  Fourier 
partial  sum  5io,o  associated  with  /  is  plotted  in  Figure  10.  To  approximate  /  by  one  of 
the  functions  Sn,m,  one  possibility  is  to  use  the  formulas  above  “blindly”  and  construct 
an  approximation,  say  59,1,  in  a  straightforward  manner.  This  approximation  is  also 
plotted  in  Figure  10.  As  the  figure  clearly  illustrates,  neither  Sio,i  nor  89,1  is  a  “good” 
approximation  to  /.  However,  we  now  observe  that  we  can  decompose  82N+2,o  into  the 
sum  of  two  other  partial  sums,  i.e. 

S2N+2.0  =  +  5)^+1, o>  where 

5^S+i,o  =  sin(x)  -  ^  sin(3x)  -I-  ^  sin(5x)  -  ..., 

■5^11.0  =  I  {^sin(2x)  -  isin(4x)  -|-  ^sin(6x)  -  ...| . 

19 


We  can  now  apply  our  formulas  to  each  of  these  sums,  separately,  and  construct  two  new 
approximations  i  and  We  then  use  +  5jv\  as  a  new  approximation  to  /. 
The  approximation  +  84^}  is  also  shown  in  Figure  10  and  is  clearly  an  improvement 
over  either  5io,i  or  59,1  . 

This  example  also  serves  to  illustrate  a  possible  advantage  of  the  more  general  form 
of  Sn,m  defined  in  equation  (4).  In  particular,  the  only  approximations  that  have  been 
investigated  here  correspond  to  setting  each  Dm  =  0  in  equation  (4).  More  generally,  the 
inclusion  of  both  sine  and  cosine  terms  in  the  denominator  of  Sn,m  allows  the  possibility 
of  “shifting”  the  location  of  an  approximate  “pole”  of  Sn,m  from  either  x  =  0  or  i  =  ±ir 
to  an  arbitrary  point  interior  to  the  interval  (— x,ir).  For  example,  when  M  =  I,  the 
denominator  of  Ss,i  can  be  expressed  as  1  +  Ci  cos(x)  +  Di  sin(ar)  =  1  +  Ci  cos(x  —  6). 
Thus,  if  Cl  —  1,  for  example,  then  the  denominator  will  be  small  when  x  is  near  6 
and,  hence,  Sn,i  could  potentially  better  simulate  a  function  which  has  a  singularity 
(or  a  “near”  singularity)  at  x  =  6. 

We  also  note  that,  while  Theorems  3  and  6  establish  the  convergence  of  the  sequence 
of  approximations  {Sn,m]  ^  fixed  value  of  A/  as  iV  — »  00,  the  Gibbs  phenomena 
has  not  been  eliminated,  although  its  effects  seem  to  be  mitigated.  Thus,  the  rate  of 
convergence  of  this  sequence  in  regions  near  a  point  of  discontinuity  of  /  needs  further 
study.  For  example,  in  the  proof  of  Theorem  3,  equation  (35)  establishes  that,  for  any 
fixed  value  of  i,  with  — tt  <  x  <  x,  the  terms  are  0(1/^”+^)  as  k  00.  However,  the 
“practical”  rate  of  convergence  of  the  series  is  mitigated  somewhat,  especially  near  x  = 
0,  by  the  presence  of  the  factor  (1  —  cos(a:))  in  the  denominator  of  the  expression  for  r^. 
This  factor,  in  combination  with  the  terms  sin{kx)  in  the  numerator,  also  foreshadows 
the  Gibbs  phenomena  that  does  remain.  However,  the  fact  that  an  explicit  form  for  these 
terms  is  available  should  assist  in  the  investigation  of  ways  to  improve  convergence  in 
these  regions.  This  observation,  as  well  as  the  other  observations,  questions,  and  several 
related  issues  raised  here,  are  the  subject  of  some  current  investigations  and  the  results 
of  these  investigations  will  be  reported  elsewhere. 
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A.  Appendix.  In  this  appendix  we  present  some  explicit  expressions  for  the  quan¬ 
tities  defined  in  section  5  in  terms  of  the  coefficients  {6*}  by  equations  (46)  and 
(49)  when  /  is  an  odd  function  of  x.  Analogous  expressions  hold  when  /  is  an  even 
function  of  x,  with  each  bk  replaced  by  a*.  For  M  =  1,2,  and  3,  these  quantities  are 
given  by 


Af  =  1  :  Si  =  — f>w+i/fciv; 


M  =  2: 


Si 


bs+ibs  —  6jv-f2^jv-i 
bs+ibs-i  —  iji/ 


32  = 


bNbN-i-2  —  bji/^i . 
&Y+1  V-1  —  br^  ' 


M  =  3: 


bN+ljb^  —  bN-ibff+l)  -I-  bs+2ibN+lbN-2  ~  bs-lbs)  +  &jV+3(fe^-l  -  bNbN-2) 

det 


bs+libN+2bN-l  -  bffbs+i)  +  bN+2{b^  —  ^^+26^-2)  +  bN+3{bN+lbN-2  —  bsbN-i) 


S3 


_  bN-i-i{b%+i  —  b^bs+T)  -b  ^•f2(fe/v-t-2fcN-i  —  bshs+i)  +  bN+3{b%  —  bN^ibs-i) 

det 


det  =  bs-2{—b%^i  +  6wV-h2)  +  bpf-i{—bN~ibN+2  +  ^bsbs^t)  -f  bsi—bh). 
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Figure  1:  A  plot  of  the  function  /(x)  (solid  line),  the  approximations  S4,o(x) 
(dashed  line)  and  53,i(x)  (dotted  line)  for  Example  1.  Note  that  the  oscillations  present 
in  ^4,0  away  from  the  singularities  of  /  have  been  virtually  eliminated  in  53,1. 


✓  ^ 


•^3,1 


Figure  2:  A  plot  of  the  function  f{x)  (solid  line),  the  approximations  5*4, o(x) 
(dashed  line)  and  53,i(x)  (dotted  line)  for  Example  2.  Note  that  the  oscillations  present 
in  84,0  away  from  the  regions  of  large  curvature  of  /  have  essentially  been  eliminated 
in  83,1. 


A^  =  2 


/4^:  4. 


iV  =  3 


Figure  3:  A  plot  of  A:^(Bfc|  (circles)  for  Example  1  as  a  function  of  1/ib,  for  2  <  .t  < 
N,  with  3  <  iV  <  20.  The  0{lfk^)  decay  of  the  coefficients  {J?*}  is  clearly  illustrated, 
as  is  an  interesting  asymptotic  behavior  of  Bn-,,  for  a  fixed  value  of  q,  as  N  oo. 
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Figure  4:  A  plot  of  A;^|Afc|/|aib|  (circles)  for  Example  2  as  a  function  of  1/A;,  for 
2  <  k  <  N,  with  3  <  TV  <  20.  The  figure  illustrates  that  the  coefficients  {A*}  decay 
faster  than  the  coefficients  {ojt},  asymptotically  by  a  factor  of  fe“*,  as  predicted  by 
Theorem  2.  As  in  Example  1,  an  interesting  asymptotic  behavior  of  As-q,  as  TV  — »  oo 
for  a  fixed  value  of  o,  is  also  predicted. 


Figure  5:  A  plot  of  log  \En,m\I  log  |A^|  for  Example  1  as  a  function  of  1/  log  |7V|  for 
3  <  AT  <  50,  at  X  =  x/2,  for  Af  =  0, 1,2,  and  3.  Here  Es^{x)  =  (/(x+)  +  /(x~))/2  - 
We  note  that  an  intercept  of  2M  +  1,  as  predicted  by  Theorems  3  and  6,  is 
clearly  consistent  with  these  plots. 


Figure  6:  A  plot  of  the  function  f{x)  (solid  line),  the  approximations  5e^(x) 
(longer  dashed  line),  55,i(x)  (shorter  dashed  line),  and  54,3(x)  (dotted  line)  for  Example 
3.  Note  that  the  oscillations  present  in  8^,0  away  from  the  regions  of  large  curvature 
of  /  have  essentially  been  eliminated  in  54,2,  while  is  not  much  of  an  improvement 
over  56,0. 
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Figure  7:  A  plot  of  log  1 5*1/ log  fc  vs.  1/logfc  for  Example  1  with  M  =  1,2,3, 
and  4,  for  N  =  20.  The  increased  rate  of  decay  to  zero  of  the  coefficients  5*  with 
increasing  values  of  M  is  clearly  evident,  especially  when  compared  with  original  Fourier 
coefficients  6*,  which  decay  only  like  1/k. 


0 


X 


TT 


Figure  8:  A  plot  of  the  exa,ct  solution  ul'*’®)  (solid  line),  as  well  as  the  approxima¬ 
tions  (dashed  line)  and  (dotted  line),  for  the  simple  heat  conduction  problem 
of  section  7,  when  i  ~  0.01  (upper  set  of  curves)  and  i  —  0.05  (lower  set  of  curves).  The 
corresponding  comparison  at  t  =  0  is  the  same  as  shown  in  Figure  1  for  0  <  x  <  t. 


Figure  9:  A  comparison  of  the  approximations  Se-g,,,  for  q  =  0,1,..., 5  for  Ex¬ 
ample  1.  In  this  case,  the  figure  seems  to  suggest  that  perhaps  is  “best"  of  the 
approximations  with  N  +  M  =  6. 
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Figure  10:  For  the  2T-periodic  function  /  defined  by  f{x)  =  xx/2,  for  — jr/2  < 
X  <  x/2,  and  f{x)  =  0,  otherwise,  in  [— jt,  x],  a  plot  of  /(x)  (solid  lines)  and  the 
approximations  5io,o  (longer  dashed  line),  89,1  (shorter  dashed  line),  and  S[]i  + 
(dotted  line)  for  0  <  x  <  ir.  The  improved  quality  of  the  approximation 
over  either  5io,i  or  89,1,  is  apparent. 
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